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Abstract 

We derive a resummation formula for a fc^-dependent parton distribution 
function at threshold, where k? is a parton transverse momentum. The 
derivation requires infrared cutoffs for both longitudinal and transverse loop 
momenta as evaluating soft gluon emissions in the Collins-Soper resummation 
framework. This unified resummation exhibits suppression at large b, b being 
the conjugate variable of kr, which is similar to the kx resummation, and 
exhibits enhancement at small b, similar to the threshold resummation. 
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1. Introduction 



Recently, we have demonstrated [jl] that both the resummation for 
a parton distribution function <p(x, kr,p + ) and the threshold resummation 
for (p(x,p + ) can be performed in the Collins-Soper (CS) framework 0. The 
distribution function kT,p + ), associated with a hadron of momentum 
pfi = p+$v+^ describes the probability that a parton carries the longitudinal 
momentum xp + and the transverse momentum kx- The distribution func- 
tion (p(x,p + ), which coincides with the standard parton model, comes from 
4>(x, k T ,p + ) with kx integrated out. In the kx resummation the double log- 
arithms ln 2 (p + 6), b being the conjugate variable of k?, are organized. The 
result is Sudakov suppression, quoted as || 

b, p + ) = exp 

where the anomalous dimension will be defined later, and (p^ is the 
initial condition of the double-logarithm evolution. In the threshold resum- 
mation the double logarithms ln 2 (l/iV), N being the moment of are 
organized. The result is an enhancement |Q: 

(p(N,p + ) = exp 

with the same anomalous dimension 7^. 

We have neglected the dependence of <fi on the renormalization (or fac- 
torization) scale n, which denotes the single-logarithm evolution, and can be 
easily derived using renormalizat ion-group (RG) equations. Equation ([l|) has 
been employed to evaluate the pr spectrum of direct photon production, pt 
being the transverse momentum of the direct photon || . It was observed that 
the Sudakov exponential, after Fourier transformed to the k T space, provides 
the necessary k^ smearing effect, which resolves the discrepancy between ex- 
perimental data and the next-to-leading-order QCD (act 2 ) predictions |J. 
Equation (|2|) is suitable for the analyses of dijet, direct photon and heavy 
quark productions in kinematic end-point regions [|J. 

It has been found that in the CS framework the two different types of 
double logarithms \n 2 (p + b) and ln 2 (l/iV) are summed by choosing appropri- 
ate infrared cutoffs in the evaluation of soft gluon corrections. If transverse 
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degrees of freedom of a parton are included, 1/6 will serve as an infrared 
cutoff for soft gluon emissions. Combined with the scale p + from the hadron 
momentum, the large double logarithms In 2 (p + b) are generated. In the end- 
point region with x — > 1, we keep the longitudinal cutoff (1 — x)p + , or p + /N 
in the moment space, and integrate out the transverse degrees of freedom of 
a parton. Combined with p + , the double logarithms ln 2 (l/iV) are produced. 
Therefore, in the case with p + b ^> N, we neglect the longitudinal cutoff, 
and sum In (p + b). In the case with iV ^> p + b, we neglect the transverse 
cutoff, and sum ln 2 (l/iV). Based on these observations, it is straightforward 
to develop a unified resummation formalism by retaining the longitudinal 
and transverse cutoff's simultaneously. It will be shown in this letter that 
this unified resummation exhibits suppression at large b, similar to the fcr 
resummation, and exhibits enhancement at small b, similar to the threshold 
resummation. 



2. Formalism 

Consider a quark distribution function for a hadron in the minimal sub- 
traction scheme, 

<j>(x,k T ,p + ) =I^J -^e-^ + y- + ^- h (p\q(y-,b)^q(0)\p) , (3) 

where 7 + is a Dirac matrix, and \p) denotes the hadron with the momentum 
p. Averages over spin and color are understood. The above definition is 
given in the axial gauge n ■ A = 0, where the gauge vector n is assumed 
to be arbitrary with n 2 ^ 0. Though this definition is gauge dependent, 
physical observables, such as hadron structure functions and cross sections, 
are gauge invariant. It has been shown that the n dependences cancel among 
the convolution factors, i.e., among parton distribution functions, final-state 
jets, and nonfactorizable soft gluon exchanges, in the factorization formula 
for a DIS structure function [Q. 

The key step in the CS technique is to obtain the derivative p + d<p/dp + 
[0]. In the axial gauge n appears in the gluon propagator, {—i/l 2 )N ia '{l), 
with 

N^(l) = J** - + y + n 2 ji±- . (4) 
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Because of the scale invariance of N^ u in n, <fi depends on p + through the 
ratio (p-n) 2 /n 2 , implying that the differential operator d/dp + can be replaced 
by d/dn a using a chain rule, 

d n 2 d 

P + TT ( f ) = v a -—(p , (5) 

dp + v ■ n dn a 

where v = (1,0,0) is a dimensionless vector along p. The operator d/dn a 
applies to N^ u , leading to 

' = v a (NT + NT) , (6) 



v ■ n dn, 
with the special vertex 



n 2 v n 



, (7) 
v ■ nn ■ I 

The momentum Z M (/") is contracted with the vertex the differentiated 
gluon attaches, which is then replaced by the special vertex in Eq. (|7|). For 
each type of vertices, there exists a Ward identity, which relates the dia- 
gram with the contraction of P (l v ) to the difference of two diagrams [[/J. A 
pair cancellation occurs between the contractions with two adjacent vertices. 
Summing the diagrams containing different differentiated gluons, the special 
vertex moves to the outer end of a parton line. We arrive at the derivative, 

described by Fig. 1(a), where the square in the new function represents the 
special vertex v a . The coefficient 2 comes from the equality of the two new 
functions with the special vertex on either side of the final-state cut. 

To obtain a differential equation of from Eq. (§), we need to factorize 
the subdiagram containing the special vertex out of 0. The factorization 
holds in the leading regions of the loop momentum I that flows through the 
special vertex. The collinear region of I is not leading because of the factor 
1 / (n-l) in v a with nonvanishing n 2 . Therefore, the leading regions of I are soft 
and hard, in which the subdiagram is factorized from <p into a soft function 
K and a hard function G, respectively. The remaining part is the original 
distribution function <p. That is, (f) is expressed as the convolution of the 
functions K and G with <b. 
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The lowest-order contribution to K is extracted from Fig. 1(b), whose 
contribution is written as 



<f) s (x, k T ,p + ) = <f) sv (x, k T ,p + ) + (f) sr (x, k T ,p + ) 



with 



ig 2 C F /i t 



ig 2 C F \i t 
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N uP (l) V -^-2m5{l 2 ) 
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(9) 



(10) 



(11) 



corresponding to the virtual and real gluon emissions, respectively. C F = 4/3 
is a color factor, and 5K an additive counterterm. The ultraviolet pole in 
Eq. (|ToD is isolated using the dimensional regularization. 

To work out the It integration explicitly, we employ the the Fourier trans- 
form from the kr space to the b space. Inserting the identities 



f 1 d£5(£ -x) = l, f 1 d^6^ -x- l + / P + ) = 1 

J X J X 

into Eqs. (|It]) and QTTD, respectively, Eq.(Q) becomes 



(12) 



<fi s (x,b,p + ) 



d£ 



Kill — 



x \ p^ 
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(13) 



with 



K = ig 2 C F \i t 



d 4 ~n 



+2m5{V)5 1 



X I 



I 2 



5K5 1 



x 



(14) 



To obtain the above expression, we have adopted the approximation 



5 1-% ' 
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which is appropriate in the threshold region with x — > 1. The exponential 
exp(il<r • b) in the second term arises form the Fourier transform of the real 
gluon contribution. 

To work out the £ integration explicitly, we further employ the Mellin 
transform from the momentum fraction (x) space to the moment (N) space: 

<j) s (N,b,p + ) = [ dxx N ~ 1 (fi s (x, b,p + ) , 
Jo 

= ^(p + /(W,&/i,a s (/i))0(iV,6,p+) . (16) 

with 

K(p + /(Nfx), bfi, a,(ji)) = ^ dzz N - l K((l - z)p + /fi, 6/x, a,(ji)) . (17) 

The convolutions between K and 4> in l + and It are then completely simplified 
into a multiplication under the Mellin and Fourier transforms, respectively. 
Equation flT7p is evaluated in the Appendix, and the result is 



K(p + /(Nfi),bn,a s (n)) = ^-C F 

7T 



bfi °l N J 



(18) 



K being the modified Bessel function. The gauge factor v = <J (v ■ n) 2 /\n 2 \ 
confirms our argument that <ft depends on p + via the ratio (p ■ n) 2 /n 2 . It 
is easy to examine the large p + b and iV limits of the above expression. For 
p + b ^> N, we have Kq — > and 

a 1 

C F hi f-, (19) 

7T b\l 

which is exactly the soft function with the characteristic scale 1/6 appearing 
in the kx resummation ]l|, For N ^> p + b, we have Kq m — hi{yp + b/N) 
and 

K^^C F ]*?£, (20) 

which is the soft function with the characteristic scale p + /N for the threshold 
resummation [J]. Hence, Eq. fllBD is indeed appropriate for the unification of 
the kr and threshold resummations. 
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The lowest-order contribution to G from Fig. 1(c) is given by 



MN,b,p + ) = G(p + /fi,a s (fi))<P(N,b,p^ 
in the b and N spaces, with 



G = -ig 2 C F fx e 



(2tt 



,4-e 



/ 



(p - : 



r + 



-5G 



(21) 



(22) 



in the limit x — > 1, where <5G is an additive counterterm. The second term, 
whose sign is opposite to that of 4> sv , is a soft subtraction. This term avoids 
double counting, and ensures a hard momentum flow in G. A straightforward 
calculation gives 



G{p + /fi,a s (fj)) 



7T 



(23) 



which is the same as the hard functions with the characteristic scale p" 1 
appearing in the and threshold resummations [l], || . 
Using (j) = (p s + cj) h , Eq. (H) becomes 



K(p + j (N/jl), bp, a s (//)) + G(p + //i, a s (ii)) 



x<fi(N, b,p^ 



(24) 



The functions and G possess ultraviolet divergences individually as indi- 
cated by their counterterms. These divergences, both from the virtual gluon 
contribution (p sv , cancel each other, such that the sum K + G is RG invari- 
ant. The single logarithms contained in K and G are organized by the RG 
equations 

H^K = - lK = -fi-fj-G . (25) 

The anomalous dimension of K, X K = fj,dSK/dfj,, is given, up to two loops, 
by! 



a 



Ik = —C F + — C 



TV 



a. 



7T 



V 36 12 18 ; 



(26) 



with rif the number of quark flavors, and Ca = 3 a color factor. 

As solving Eq. (j25D , we allow the variable /i evolves from the characteristic 
scale of K to the scale of G. We discuss the cases for p + b ^> N and for 
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iV ^> p + b first, which will help the derivation of the unified resummation. 
The solution of K + G is written as 

K(p + 1 (Nfi), b/i, a s ((j,)) + G(p + /fi, a s (/j,)) 



1/6 jl 



d\i 



for p + b > iV , 
for N > p + b 



"- lK {a s {n)) , for iV»p+&, (27) 

lp+/N fj, 

where the initial conditions of K and G of the RG evolution have been ne- 
glected, since they are irrelevant to the double-logarithm summation. Equa- 
tion ( p7|) indicates that the distribution function (p(N,b,p + ) involves ln(p + 6) 
and ln(l/iV) in the p + b — > oo and N oo limits, respectively, as stated 
before. This can be easily understood by ignoring the variation of 7^, and 
performing the \x integration directly. 

Inserting Eq. fl2"T| ) into (|2^), we obtain the solutions 

fP + dp /-p dfi 



<f>(N, i>,p H 
0(iV, 6,p H 



exp 



exp 



-2 



1/6 p Ji/b /i 
rP + /N dp fP + dfi 



-7tf(a a (//)) 



Jp+ p Jp jj, 



(28) 
(29) 



for p + b 3> iV and iV ^> p + b, respectively. Note that we have replaced the 
derivative p + d/dp + by (p + /N)d/d(p + /N) as deriving Eq. (|29|) 0, since we 
intend to resum ln(l/7V). Obviously, Eq. ( |2"g| ) is identical to Eq. (|I]) with 
x — > 1, and Eq. ( p9|) is equivalent to Eq. (0) up to corrections suppressed by 
a power 1/iV |L|. 

Equation (0) implies the characteristic scale of i\" for the unified resum- 
mation, 



1 



exp 



p + b 



where the gauge factor 1v has been suppressed. Hinted by Eqs. 
we rewrite Eqs. (B3) and (EW) in terms of the above scale as 



<f>(N,b,P + ) 
<f>(N,b,p+) 



exp 



exp 



dp fP dfj. 



exp[-2fo(p+&/W)]/6 P Jl/b fi 

dp fP d\i 



exp[-K (p+b/N)]/b p Jp+ /i 



7K-(a s (/i)) 



(30) 
and 

(31) 
(32) 
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At last, to unify the above expressions, we replace the lower bounds of \i by 

^exp[-K (p + b)} , (33) 
which is motivated by Eq. (BO). The unified resummation is then given by 



4>(N, b,p + ) = exp 



dp fP dn 

exp[-K (p+b/N)]/b p Jex P [-K (p+b)]/b fi 



which is appropriate for arbitrary p + b and N. It is easy to justify that 
Eq. ( 151 ) approaches the kx resummation in Eq. fl28p as b — > oo, and ap- 
proaches the threshold resummation in Eq. (B^) as b — > 0. 



3. Discussion 



Comparing Eq. fl34|) with ( p28|) and fl29|) , we have the following remark. 
In the /ct resummation it is the logarithms of the large scale p + from G that 
are organized, and the result is a suppression. In the threshold resummation 
it is the logarithms of the small scale p + /N from K, as shown in Eq. (|20|), 
that are organized. The result is an enhancement. The opposite effects 
of these two resummations are attributed to the opposite directions of the 
double-logarithm evolution: from 1/b to the large p + in the former case, 
and from p + to the small p + /N in the latter case. The unified resummation 
for a /^-dependent parton distribution function at threshold exhibits both 
behaviors: it is a suppression in the large b region (p + b 3> N), and turns into 
an enhancement in the small b region (N ^> p + b). That is, Eq. (|34] ) displays 
the opposite effects of the kx and threshold resummations at different b. 

The behavior of the unified resummation can be explained as follows. For 
an intermediate x, virtual and real soft gluon corrections cancel exactly in 
the small b region, since they have almost equal phase space. Hence, there 
are only single collinear logarithms, namely, no double logarithms. In this 
case the Sudakov exponential approaches unity as b < l/p + ||, indicating 
the soft cancellation stated above. However, at threshold (x — > 1), real 
gluon emissions still do not have sufficient phase space even as b — > 0, and 
soft virtual corrections are not cancelled exactly. In this case the double 
logarithms ln 2 (l/iV) persist and become dominant. Sudakov suppression 
then transits into an enhancement, instead of unity, smoothly as b decreases. 
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The result obtained in this work will be useful for the study of, say, the 
production of large pr direct photons or jets. Phenomenological applications 
of Eq. ( P^D will be published elsewhere. 

This work was supported by the National Science Council of Republic of 
China under Grant No. NSC-88-2112-M-006-013. 



Appendix 



In this Appendix we present the detail of deriving Eq. (|18|) . A simple 
investigation of Eq. (Q) shows that the first term g v p in N u/3 gives a vanishing 
contribution because of v 2 = 0, and the contribution from the second term 
—n v lp/n-l cancels the contribution from the fourth term n 2 /^/ (n-l) 2 . Hence, 
we concentrate only on the third term —npl v /n ■ Z, which leads Eq. (|T7j ) to 



K 



-ig 2 C F ^ / dzz N - 1 
Jo 



n 



(2vr) 4 - e (n ■ iy 



6(1 



I 2 



+ 27TiS{l 2 )S 1 



P 



- 5K6(1 - z) 



(35) 



Performing the integrations over I and Z + , we have 



K 



—g Cp\i e \ dzz 



JV-l 



10 

2n 2 (l 



d 2 ~H 



(2< 



3-e 



5(1 



1/rp 



z)p 



+2 



[n+l 2 T + 2n~(l -z) 2 p+ 2 } 2 



5K5{1 



The first term in the integral has been evaluated in 
with its ultraviolet pole subtracted by 5K is quoted as 



(36) 



and the result 



^ln^ 

7T jJL 



(37) 



where constants of order unity have been neglected. The infrared regulator a, 
introduced by the replacement of the denominator l\ by l^ + a 2 , will approach 
zero at last. The second term, after performing the integration over It, gives 



2^C F vp + b jT 1 dzz N - l K 1 (2^/(l 



z) 2 u 2 p +2 + a 2 b 



(38) 
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where K\ is the modified Bessel function, and the infrared regulator a for 
the divergence at z — > 1 is also introduced by the same replacement. 
We adopt the identity 



£ dzz N - l K l (2y/(l - zfu 2 p+ 2 + a 2 b 

= f 1 dz(z N - 1 - l)Kx(2(l - z)vp + b) 
Jo 

+ J dzK x (2yJ(l - z) 2 v 2 p+ 2 + a?b 



(39) 



where a in the first term has been dropped, since the ^-integral is infrared 
finite. To perform the integration, we employ the relation 



dz{z N - 1 - 1)^(2(1 - z)vp + b) 



l-l/N 



dzK 1 (2(l- z)vp + b) , (40) 



which is is valid up to corrections suppressed by 1/JV. It is then trivial to 
show that Eq. ( |3~8|) leads to 



a, 



c f 



7T 



\n--K 

ab 



'2vp+b s 
N 



(41) 



At last, combining Eqs. (|37|) and (PD, we obtain Eq. (p!8[). 
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Figure Captions 

Fig. 1. (a) The derivative p + d(j)/dp + in the axial gauge, (b) The 0(a s ) 
function K. (c) The 0(a s ) function G. 
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